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We review different models used for reactions involved in nuclear astrophysics. The 
reaction rate is defined for resonant as well as for non-resonant processes. For low- 
density nuclei, we describe the DWBA method, the potential model, the i?-matrix 
method, and microscopic cluster models. The statistical model is developed for 
high-level densities. Details of calculations in the low- and high-density regimes are 
illustrated with new results concerning transfer reactions and level densities. 

O 

P ! 1 Introduction 

h 

Nuclear reactions are important for studying energy generation and nucleosyn- 
thesis in astrophysical processes [1,2]. Especially in explosive burning scenar- 
ios, large reaction networks include mainly nuclei far from stability which have 
not been studied experimentally so far, and perhaps will never be accessible 
in the laboratory. This poses a challenge to both experimentalists and espe- 
cially theorists, which have to predict reaction cross sections of more or less 
unknown nuclei. Another challenge is found in the comparatively low interac- 
tion energies which pose problems for predictions as well as for measurements 
already close to and at stability. Here we provide and overview of the current 
status of theory and details of the remaining challenges, illustrated by several 
examples of new calculations. 

Theoretical models can be roughly classified in three categories [3]: 
(i) Models involving adjustable parameters, such as the i?-matrix [4] or the 
X-matrix [5] methods; parameters are fitted to the available experimental 
data and the cross sections are extrapolated down to astrophysical energies. 
These fitting procedures of course require the knowledge of data, which are 
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sometimes too scarce for a reliable extrapolation. 

(ii) "Ab initio" models, where the cross sections are determined from the 
wave functions of the system. The potential model [6], the Distorted Wave 
Born Approximation (DWBA) [7], and microscopic models [8,9] are, in princi- 
ple, independent of experimental data. More realistically, these models depend 
on some physical parameters, such as a nucleus-nucleus or a nucleon-nucleon 
interaction which can be reasonably determined from experiment only. The 
microscopic Generator Coordinate Method (GCM) provides a "basic" descrip- 
tion of a A-nucleon system, since the whole information is obtained from a 
nucleon-nucleon interaction. Since this interaction is nearly the same for all 
light nuclei, the predictive power of the GCM is important. 
(Hi) Models (i) and (ii) can be used for low level-density nuclei only. This 
condition is fulfilled in most of the reactions involving light nuclei (A < 20). 
However when the level density near threshold is large (i.e. more than a few 
levels per MeV), statistical models, using averaged optical transmission coef- 
ficients, are more suitable (see Sec. 4). 

In Sec. 2, we start with useful definitions commonly used in nuclear astro- 
physics, such as of the reaction rates which are the actual quantities used in 
reaction networks, and their relation to cross sections and the relevant en- 
ergy range. We then discuss different approaches and reaction mechanisms for 
systems with low (Sec. 3) and high level densities (Sec. 4). In nuclear astro- 
physics, the relevant projectile energy range often covers both the low and 
high density regimes. The difficulties concerning the prediction of reaction 
rates bridging the two regimes are briefly addressed in Sec. 4.3. Concluding 
remarks are presented in Sec. 5. 



2 Reaction Rates 



Let us consider a reaction between two nuclei with masses Aiin^ and A 2 rriN 
and charges Z\e and Z 2 t (we express here the masses in units of the nuclear 
mass m^). The theory of stellar evolution involves the reaction rate at tem- 
perature T, defined as [10,1,2] 



N A <av> = N A ( tWvT] [v(E)Eexp(-E/k B T)dE, (1) 

\TT/j,m N (k B T) s J J 

where we assume that the star can be considered as a perfect gas following 
the Maxwell-Boltzmann distribution. In Eq. (1), v is the relative velocity, 
Na is Avogadro's number, ji is the dimensionless reduced mass, and k B is 
Boltzmann's constant. At sub-coulomb energies the cross section varies as 
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a(E) ~exp(-2irri)/E, (2) 
where 77 = Z 1 Z 2 e 2 /hv is the Sommerfeld parameter. 

Using (2), the integrand of (1) can be approximated by a Gaussian shape 
[10,1] with a maximum at the Gamow peak. The energy and width are given 
by 



e 2 



7T—Z 1 Z 2 k B T(fxm N c 2 /2) 1 / 2 



2/3 

« 0.122 /2 1/3 (Z 1 Z 2 T 9 ) 2/3 MeV, (3) 



A£ = 4(E k B T/3)^ « 0.237 {Z\Z\\i) x l % T 9 5/6 MeV, (4) 

where T 9 is the temperature expressed in 10 9 i^. The Gamow energy defines the 
energy range where the cross section needs to be known to derive the reaction 
rate. In most cases, this energy is much lower than the Coulomb barrier which 
means that the cross sections drop to very low values. To compensate the 
fast energy dependence of the cross section, nuclear astrophysicists use the 
^-factor defined as 



S(E) = a(E) E exp(27r?y), (5) 

which is mainly sensitive to the nuclear contribution of the cross section. For 
non-resonant reactions, the energy dependence of the S factor is smooth. 

Rigorously the reaction rate should be calculated numerically by using experi- 
mental or theoretical cross sections. However, the analytical approach provides 
a more intuitive understanding of the physics, and is still widely used. Let us 
start with non-resonant reactions, where the S'-factor weakly depends on en- 
ergy. In this case, the integral (1) can be replaced by an accurate analytical 
approximation. A Taylor expansion near E provides 



exp(-27r?7 - E/k B T) w exp(-3E /k B T) exp [ - ( ) • (6) 

Assuming a linear variation of S(E) in the Gamow peak [1], the reaction rate 
is then given by 



Ar / 32£ \* / 3£ \ A 5k B T" 
N A <av>^N A [ - exp -— £ 1 + 



\3fj,mN(k B T) 3 J \ k B T J \ 36E J 
xS(E + h B T), (7) 



which presents a fast variation with temperature. 
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The formalism developed here is valid for any form of the S-factor. The an- 
alytical procedure can be extended further by assuming a quadratic form of 
the S'-factor 



S(E)^S + S' E + ^E 2 , (8) 



which is used in some astrophysics tables [10]. We have 



, T ( 32 ^o V ( 3£ 
Na < cry >Ri N A exp 



X 



\3fj,m N (k B T) 3 J \ k B T 

CjL T -1 OCT i on 

Soil + ||r) + S' (E + ^k B T) + -S'*E G {E + f Q k B T) 



(9) 



This yields the well known T 1//3 expansion of the reaction rate, up to T 5//3 [10]; 
Eq.(9) provides the rate from the S'-factor properties at zero energy. 

For resonant reactions, the S-factor is assumed to have a Breit-Wigner form 
(for the Breit-Wigner expression for cross sections see Eq. 28). The general 
definition (1) is of course still valid. However, one has to account for the fast 
variation of S(E) near the resonance energy. Since a numerical approach is 
difficult for narrow resonances, we present an analytical method, widely used 
in nuclear astrophysics. 

A careful analysis of integrand (1) shows that it always presents two maxima: 
at the resonance energy, and at the Gamow energy. The peak at the resonance 
energy does not depend on temperature. The second peak, corresponding to 
the Gamow energy, moves according to the temperature. From these consid- 
erations, and except in the temperature range where both peaks overlap, the 
resonant reaction rate can be split in two terms 

N A < av >^ N A <crv > R +N A < av > T , (10) 

where N A < crv > B corresponds to the maximum at E — Er. For a narrow 
resonance, we have 



N A < av > R = N A [ — H u'y ex.p(—E R / k B T) , (11) 

\Hm N k B T j 

where the resonance strength wy is defined by 



2j + i rir 2 
U1 ~ (2/ 1 + i)(2/ 2 + i) fT+ry (12) 
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J being the spin of the resonance, and (Y\, T 2 ) the widths in the entrance and 
exit channels. In capture reactions, the 7 width is in general much lower than 
the particle width. The resonance strength is then proportional to the lower 
partial width. 

The second maximum of integrand (1) yields the so-called "tail resonance" 
term Na < o~v >r- Its analytical expression is identical to the non-resonant 
rate (7) with a Breit-Wigner expression for S(E). 



3 Systems with low level densities 

3. 1 Introduction 

Although it is in general not trivial to draw a distinct line between the direct 
and the compound reaction mechanism (as, e.g., in the case of processes in- 
volving only few interaction steps of a few nucleons without participation of 
the bulk of nucleons), direct reactions are usually assumed to be fast, one-step 
processes with few degrees of freedom in which the captured particle directly 
enters the final state without sharing any energy with other nucleons [11]. It 
has long been known that direct reactions mainly occur at high projectile ener- 
gies where the formation of a compound nucleus is suppressed due to the slow 
reaction time-scale [11,12]. Such high projectile energies (-E pro j > 20 MeV) 
are not important in the effective energy range relevant for hydrostatic or 
explosive astrophysical burning. However, direct reactions can also dominate 
the cross section at very low energies in systems with low level densities. This 
occurs for light nuclei and for nuclei with low particle separation energies. For 
example, very neutron-rich nuclei far off stability have decreasing neutron sep- 
aration energies which implies that low level densities are encountered in the 
projectile-target system at r-process energies. Similar considerations hold for 
proton-rich nuclei close to the proton dripline encountered in the rp-process. 
A special case concerns nuclei with magic nucleon numbers. Due to their wide 
level spacing, direct reactions may become important for intermediate and 
heavy nuclei already close to or at the line of stability (see also Sec. 4.2). On 
the other hand, reactions in light nuclei are always direct because the reaction 
systems contain only few nucleons. Due to this fact, reactions in solar burning 
have been known to proceed via the direct mechanisms since decades. More 
recently, a number of studies have analyzed reactions in astrophysical burning 
(for both capture and transfer reactions) and also investigated direct inter- 
actions in neutron capture on intermediate and heavy nuclei far off stability 
occurring in the r-process [13,14,15,16,17]. 

In fact, the distinction between direct and compound reactions might be mis- 
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leading, especially in light nuclei. For example, it is often assumed that the 
term "direct" is synonymous to "non-resonant". This is incorrect since po- 
tential resonances can still be included in the models described below (e.g. 
[13]). Only very narrow single-particle resonances are not reproduced in mod- 
els using effective optical potentials. Therefore we prefer to address reactions 
in low level density systems and such in systems with high level densities, and 
to present the relevant models for cross section predictions in this context. 
Nevertheless, for simplicity we will still occasionally use the term "direct" in 
the following discussion. However, it is to be understood as described above. 

The discussion of reactions in low level-density systems is started with dis- 
torted wave methods for transfer and capture reactions, belonging to the group 
of potential models which describe the dynamics of a system by solving the 
Schrodinger equation with an effective potential. The potential is of optical- 
potential type (see also Sec. 4.1) but other nucleus-nucleus potentials are used 
in different models due to the nature of the model and varying approximations 
of coupling to other channels [11]. 

3.2 The DWBA method for transfer reactions 

The Distorted Wave Born Approximation (DWBA) starts from the premise 
that elastic scattering is dominant and has to be treated fully, while non-elastic 
events can be treated by perturbation theory. Although DWBA is a first-order 
theory, the way it is usually applied is not simply first-order. That is because 
optical potentials fitted to elastic scattering data may include higher-order 
effects implicitly. Therefore different potentials are needed for higher-order 
methods, such as coupled-channels calculations, than those used in DWBA in 
order to reproduce the same elastic data. 

The DWBA method can be applied to transfer reactions 

a(— b + x)+A—>b + B(=A + x) (13) 

and assumes that particle x goes from the projectile a to the target A [11]. 
The cross section for reaction (13) is obtained from the matrix elements 

t dwba = J J dpadp ^ {pp) < ^ B \AV\^ a ^ A > x +(pJ, (14) 

where the distorted waves Xa(Pa) an d Xp(Pp)i (Pa anc ^ Pp are the relative 
coordinates) corresponding to the relative motion in the entrance and exit 
channels, respectively, are generated by optical potentials U a and Up. The 
residual interaction is defined in two different ways 
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AV = V xA + V bA -U a (prior) 
= V bx + V bA -U p (post), 



(15) 



which correspond to "prior" and "post" definitions, respectively; they provide 
identical values for T DWBA . The main problem of the method is that the 
potentials are usually poorly known. In general, a good approximation is to 
neglect V bA — U a or V bA — Up. The matrix element then contains distorted 
scattering wave functions Xa, XPi an d the radial bound state wave functions 
of the transferred cluster. All of these wave functions are usually numerically 
computed with optical potentials as shown in Sec. 4.1. 

Since more realistic descriptions of nucleus a(B) should involve other configu- 
rations than b + x (A + x), spectroscopic factors are introduced (<S a and Sb)- 
The DWBA cross section is therefore linked to the experimental cross section 
through 

Ocxp = <S a <SB o"dwba- (16) 

A recent work [18] aims at investigating the precision of the DWBA for transfer 
reactions at low energies. We have investigated the 13 C(a,n) 16 reaction by a 
microscopic model, and used the DWBA in conditions as close as possible to 
the reference calculation. The conclusion is twofold. On one hand, the DWBA 
method turns out to be very sensitive to the conditions of the calculations: 
choice of the nucleus- nucleus potentials and, to a lesser extent, of the internal 
wave functions of the colliding nuclei. This sensitivity is due to very basic 
properties, i.e. the short-range character of the DWBA matrix elements, which 
are quite sensitive to details of the wave functions. On the other hand, the 
difference between the DWBA and the reference microscopic method can be 
fairly large, and varies with angular momentum. This is most likely due to 
antisymmetrization effects which are approximately included in the DWBA 
through the choice of deep nucleus-nucleus potentials. This property should 
also occur in other systems and suggests that the DWBA method can only 
provide transfer cross sections with a non-negligible uncertainty. 

Future experiments with unstable beams will provide valuable information on 
the level structure and spectroscopic factors (also for the important neutron 
capture reactions discussed in Sec. 3.3). Usually, (d,p) reactions are used for 
this purpose as neutron capture at low energies cannot be studied on radioac- 
tive isotopes. However, an analysis of the results might be complicated because 
the experimental cross sections have to be compared to calculations which still 
bear considerable uncertainties in the optical potentials used. One has to be 
extremely cautious in interpreting such comparisons as the resulting spectro- 
scopic factors may carry large uncertainties. A promising way out might be 
the use of the adiabatic approximation instead of DWBA which arises from 
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neglecting the internal excitation energies of states relative to the bombarding 
energy [11,19]. It can be formulated in such a manner that no optical poten- 
tials are needed [20]. Some parameters have to be adjusted by including elastic 
data but no further reaction information is required. Further investigation will 
show whether this approach is practicable for reactions far off stability. 



3.3 Capture reactions in a potential model 



Capture reactions can be treated in a first-order approach involving an electro- 
magnetic operator describing the emission of 7-radiation due to the dynamics 
in the movement of the electric charges. Since the application to neutron or 
proton capture dominated by electric dipole transitions is most important in 
astrophysics, here we give the expression for El direct capture. The cross sec- 
tion of the direct radiative capture process a + A — * 7 + B, which is entirely 
electromagnetic, is given by [2,21,22]: 



16tt /£ 7 \ 3 e 2 /4 

a El = IT" 



9 \hcj h 2 k a (2/ a + 1) (27a + 1) \m. 
x J2 (2^+l)(2J a + l)max(4,^) 



z a _ Za y 



C 2 S, 



J (3 



(17) 



X 



t/3 l a 



1 2 I D 1 2 

> aj\Kip a \ 



I J a J/3 

with the radial integral 

Ripa = JT f u }( r )x a {r) r dr 



(18) 



In the above expression, the energy of the emitted photon is E 1 . The charge 
and mass of the projectile and target nucleus are Z a , m a , Z A and m^, re- 
spectively. The orbital and total angular momentum quantum numbers of the 
nuclei in the entrance and exit channels are £ a , J a , £p and Jp, respectively. 
The spin quantum number, the orbital and total angular momentum quantum 
numbers are characterized by S, L and I, respectively, with indices a, A and 
B corresponding to the projectile, target and residual nucleus, respectively. 
The symbol { . . . } is the Qj symbol. The radial distorted wave in the entrance 
channels is given by Xa- The radial bound state wave function of the final state 
is denoted by up. The spectroscopic factor and the isospin Clebsch-Gordan 
coefficient for the partition B = A + a are given by C and Se^jp, respectively. 

The required wave functions Xa, up are found by solving the Schrodinger 
equation with an appropriate optical potential. The spectroscopic factor is the 
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overlap of the wavef unctions of the target and the final nucleus and is extracted 
from experiment or a more microscopic calculation. In case of single nucleon 
capture the spectroscopic factor is simply related to the occupation probability 
V 2 of the orbital into which the nucleon is captured, e.g., S = (1 — Vj) 
for depositing a nucleon onto an even target [12]. The same spectroscopic 
factor applies to (d,p) and (d,n) reactions, respectively. For a discussion of how 
to extract spectroscopic factors for neutron or proton capture such transfer 
reactions, see the end of the previous Sec. 3.2. 

The optical potential can be derived using the folding approach (see Sec. 4.1) 
and by determining the parameter A from systematics of volume integrals [23]. 
The capture cross sections are very sensitive to excitation energy and spin of 
the final states. This poses a big problem far off stability where there is no 
experimental information available. It has been shown that neutron capture 
cross sections predicted utilizing input from different microscopic structure 
models can differ by orders of magnitude [24]. Large jumps in the cross sec- 
tion when going along an isotopic chain were also found, due to the fact that a 
weakly bound state (which nevertheless dominates the cross section) can be- 
come unbound in the subsequent isotope and become unavailable for capture. 

3.4 R-matrix 

In the .R-matrix method, the energy dependence of the cross sections is ob- 
tained from Coulomb functions, as expected from the Schrodinger equation. 
The goal of the .R-matrix method [4] is to parameterize some experimentally 
known quantities, such as cross sections or phase shifts, with a small number 
of parameters, which are then used to extrapolate the cross section down to 
astrophysical energies. 

The .R-matrix framework assumes that the space is divided into two regions: 
the internal region (with radius a), where the nuclear force is important, and 
the external region, where the interaction between the nuclei is governed by 
the Coulomb force only. Although the .R-matrix parameters do depend on the 
channel radius a, the sensitivity of the cross section with respect to its choice 
is quite weak. The physics of the internal region is parameterized by a number 
iV of poles, which are characterized by energy E\ and reduced widths 7a«. In 
a multichannel problem, the .R-matrix at energy E is defined as 

N 

R aP {E) = £ (19) 

A=l A 

which must be given for each partial wave J. Indices a and (5 refer to the 
channels. 
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Fig. 1. i?-matrix fit of the 14 N(p,7) 15 reaction for capture to the 15 ground-state. 
Open circles are the data of Schroder et al. [29] and full circles represent the LUNA 
data [28]. 

The i?-matrix method can be applied to transfer as well as to capture reac- 
tions. It is usually used to investigate resonant reactions but is also suited 
to describe non- resonant processes [25]. In the latter case, the non-resonant 
behavior is simulated by a high-energy pole, referred to as the background 
contribution, which makes the R- matrix nearly energy independent. The pole 
properties are associated with the physical energy and width of resonances, 
but not strictly equal. This is known as the difference between "formal" pa- 
rameters and "observed" parameters, deduced from experiment. In a general 
case, involving more than one pole, the link between those two sets is not 
straightforward; recent works [26,27] have established a general formulation 
to deal with this problem. 

Figure 1 shows a recent application to the 14 N(p,7) 15 reaction (ground-state 
contribution) with the data obtained by the LUNA collaboration [28] . Within 
the El multipolarity, l/2 + and 3/2+ resonances are taken into account. The 
ground-state contribution represents about 15% of the total S factor. At zero 
energy, we have S(0) = 1.7 ± 0.2 keV-b, in agreement with previous results. 

3.5 Microscopic models 

Microscopic models are based on basic principles of quantum mechanics, such 
as the treatment of all nucleons, with exact antisymmetrization of the wave 
functions. The hamiltonian of an A-nucleon system is 
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i=l i<j=l 



(20) 



where Tj is the kinetic energy and a nucleon-nucleon interaction [8,9]. 

The Schrodinger equation associated with this hamiltonian can not be solved 
exactly when A > 3. The Quantum Monte Carlo method [30] represents a 
significant breakthrough in this direction, but is currently limited to A = 10. 
In addition its application to continuum states is not feasible for the moment 
(it has been applied to the d(a,7) 6 Li reaction [31], but the a + d relative 
motion is described by a nucleus- nucleus potential). 

In cluster models, it is assumed that the nucleons are grouped in clusters. We 
present here the application to two-cluster systems. This provides a natural 
extension of the formalism to scattering states. The internal wave functions of 
the clusters are denoted as 0[ i7W (&) , where Ii and 7Tj are the spin and parity 
of cluster i, and $,i represents a set of its internal coordinates. One defines the 
channel function as 



JM 



(21) 



where different quantum numbers show up: the channel spin I, the relative 
angular momentum £, the total spin J and the total parity n = 7r 1 7r 2 (— Y . 



The total wave function is written as 



= Y.Agiup) ^rovad), (22) 

all 

which corresponds to the Resonating Group (RGM) definition. Index a corre- 
sponds to different two-cluster arrangements, and A is the antisymetrization 
operator. In most applications, the internal cluster wave functions are de- 
fined in the shell model. Accordingly, the nucleon-nucleon interaction must be 
adapted to this choice, which leads to effective forces, such as the Volkov [32] 
or the Minnesota [33] interactions. The relative wave function g(p) is to be de- 
termined from the Schrodinger equation. In recent applications, this function 
is expanded over gaussian functions [8,34], which corresponds to the Genera- 
tor Coordinate Method (GCM). The GCM is equivalent to the RGM, but is 
better adapted to numerical calculations, as it makes uses of projected Slater 
Determinants. 

The main advantage of cluster models with respect to other microscopic the- 
ories is its ability to deal with reactions, as well as with nuclear spectroscopy. 
The first applications were done for reactions involving light nuclei, such as 
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Fig. 2. 17 F(p,7) 18 Ne astrophysical 5-factor with the contribution of different mul- 
tipolarities. 

d, 3 He or a particles. More recently, much work has been devoted to the im- 
provement of the internal wave functions: multicluster descriptions [35,36], 
large-basis shell model extensions [37], or monopolar distortion [38]. Micro- 
scopic cluster models have a wide range of applications, both in low-energy 
reactions and in the spectroscopy of light nuclei. The main limitation arises 
from the number of channels included in the wave function, which reduces the 
validity of the model at low energies. Also high-level densities require many 
channels in the wave functions. Nuclear astrophysics is probably one of the 
best candidates for applications of microscopic cluster models. The low en- 
ergies and the low densities involved in most nuclei make the conditions of 
application quite valid. 

We illustrate cluster models with a recent application to the 17 F(p,7) 18 Ne 
reaction [39]. The knowledge of the 17 F(p,7) 18 Ne reaction rate is important 
for the understanding of the energy production and of the nucleosynthesis 
in novae [40]. Until now, a direct measurement of the cross section down to 
astrophysical energies (typically 0.3 MeV) has not been performed. It is now 
well established that the 3f level dominates the rate at high temperatures since 
it is the only s— wave resonance at low energies. The gamma width remains 
experimentally unknown and is estimated by shell-model calculations. At low 
temperatures, the non-resonant contribution is dominant but is evaluated with 
simple models [41] which need confirmation from more elaborate theories. 

The GCM wave functions are defined as a combination of 17 F + p and «+ 14 
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channel functions. The 17 F and the 14 nuclei are described by s,p and sd 
(for 17 F) harmonic oscillator shells. This framework allows the description of 
the 17 F(5/2+ 1/2+ 3/2+)+p and 14 O(0+2+) +a channels. 

The proton width of the 3 + resonance (21 keV) is in fair agreement with the 
experimental value (18 ± 3 keV) [42]. For the gamma width, we find 33 meV; 
a recent evaluation of the reaction rate [42] uses T 7 = 29 ± 19 meV. The 
astrophysical S'-factor is calculated for the El, E2, and Ml multipolarities. 
In each case, the four 0+,2 + , 2 + and 4 + bound states are considered. The J 
values for scattering states are limited to 4, which corresponds to £ = 2. The 
total and partial ^-factors are displayed in Fig. 2. 



4 Systems with high level densities 

4-1 Optical model 

Due to the complexity of the nucleon-nucleon (NN) interaction, one often 
resorts to working with effective interactions instead of solving microscopic 
models based on NN potentials (but see Sec. 3.5). Widely used in calculating 
different reaction mechanisms is the optical model. In that model, the compli- 
cated many-body problem posed by the interaction of two nuclei is replaced by 
the much simpler problem of two particles interacting through an effective po- 
tential, the so-called optical potential [11,12,43]. Such an approach is usually 
feasible only with few contributing channels. Always included is elastic scat- 
tering. That is why optical potentials can be derived from elastic scattering 
data. 

In general, the optical potential is complex. The real part describes the refrac- 
tion of the incident waves by the nucleus (elastic scattering), the imaginary 
part their absorption, hence the term "optical" model. Thus, the imaginary 
part of the resulting wave function contains all non-elastic processes remov- 
ing flux from the incident channel. In general, the potential is expected to be 
energy-dependent for obvious reasons. Most applications make use of local po- 
tentials although, in principle, also the effective potential should be non-local. 
An additional slight energy-dependence is introduced by this neglection [11]. 

The required potential parameters can be derived either by fits to experimental 
elastic data or from microscopic considerations. An example of the latter is the 
folding approach [44] which is often employed to determine the real potential. 
In this approach an effective NN interaction is folded with projectile and target 
nucleon density distributions. The latter have to be derived from electron 
scattering data or further microscopic considerations. Another example for 
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a more microscopic approach is the use of the local density approximation 
[45,46]. The advantage of this model is that also the imaginary part can be 
determined. Again, certain parameters have to be fitted to experimental data 
and density distributions have to be known. 

The standard definition of the optical model stems from the description of elas- 
tic and inelastic scattering. The time-independent radial Schrodinger equation 
is numerically solved with an optical potential which provides a mean inter- 
action potential, averaging over individual NN interactions. In consequence, 
single-particle resonances cannot be described in such a model. However, res- 
onances stemming from potential scattering can still be found. Elementary 
scattering theory yields expressions for the elastic cross section and the re- 
action cross section. The latter includes all reactions and inelastic processes 
which cause loss of flux from the elastic channel. With the diagonal element 
S aa of the S-matrix S (sometimes also called scattering matrix or collision 
matrix), the reaction cross section for spinless particles is then given by [43] 

< Q = |rE(2£ + l)(l-m 2 ) • (23) 

This can be generalized to other channels f3, not just the elastic one. The 
elements of the S-matrix are complex, in general, and related to the scattering 
amplitude / of the outgoing wave function 

r /3 = 4^ (2L + 1)( ^ /3 " 1)PL(cos ^ ) ' (24) 

L 

which is nothing else than the transition amplitude T Q/3 = — (27rft)//i /3 /' a/3 , 
in turn, with (j, being the reduced mass. As mentioned above, the imaginary 
part of the optical potential gives rise to an absorption term in the solution 
of the Schrodinger equation, thus removing flux from the considered channels. 
Therefore, the matrix element S aa is also related to the transmission coefficient 

2 ) = 1 - \St a \ 2 (25) 

which describes the absorption of the projectile by the nucleus. Important for 
practical application is that the phase shifts 5l can be derived from elastic 
data. 

In the preceding Sec. 3.1 we have already discussed several models making 
use of potentials of the optical type. In the following, the statistical model of 
nuclear reactions is described which also applies optical potentials to calculate 
the formation and decay of the compound nucleus. 
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4-2 Statistical model 



It is well established that the number of states in a nucleus increases in prin- 
ciple exponentially with excitation energy. Thus, with increasing projectile 
energy, different resonance regimes are covered. Starting without or only few 
resonances at low energy and more, but still well separated, resonances at in- 
termediate energies, one will eventually encounter a region where the average 
resonance width (r) becomes larger than the average level spacing D = l/p. 
This is the domain of the compound nucleus reaction which can be calculated 
in a statistical model only using average resonance properties, also called the 
Hauser-Feshbach approach. In astrophysical applications, the compound reac- 
tion mechanism is especially important for reactions with charged projectiles 
as the relevant energy range (given by the Gamow peak, eq. 4) is shifted to 
higher energies as compared to neutron-induced reactions. Nevertheless, also 
for reactions involving neutrons the statistical model can be applied when the 
reaction Q-value or the level density is sufficiently high. This is true for most 
intermediate and heavy nuclei close to stability with non-magic neutron num- 
bers. The relevant quantity is the number of available levels in the Gamow 
peak. An estimate of the applicability range of the statistical model is given 



The statistical model is based on the Bohr independence hypothesis [48]. It 
states that the projectile forms a compound system with the target and shares 
its energy among all of the nucleons. This implicitly requires long reaction 
time-scales as the compound nucleus has to live long enough to establish com- 
plete statistical equilibrium. Compared to the direct mechanism (Sec. 3.1) it 
is about 5 orders of magnitude longer and includes many degrees of freedom. 
Finally, the compound nucleus can decay by emitting photons or particles. Ac- 
cording to the Bohr hypothesis the final decay is independent of the formation 
process, except for the basic conservation laws of energy, angular momentum, 
parity, and nucleon number. Thus, the cross section can be factorized into two 
terms 



the formation cross section <7f orm and a branching ratio describing the prob- 
ability for decay to the observed channel (3. An early implementation of this 
was the Weisskopf-Ewing theory [49]. Today, the Hauser-Feshbach approach 
[50] is widely used which also incorporates angular momentum conservation. 
Modern theories also treat pre-equilibrium emission of particles in different 
ways. However, for the energies encountered in astrophysical scenarios, this 
additional mechanism is usually negligible. 



in [47]. 



^CN — °form°dcc — Oft 




(26) 



orm 



(r tot > 
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In the Hauser-Feshbach theory, the formation cross section is calculated as 
in the optical model (Sec. 4.1) but using averaged transmission coefficients 
(T) which can again be computed from a Schrodinger equation and using 
appropriate optical potentials. The average widths occurring in the branching 
ratio &dec are related to average transmission coefficients via (T) = 2np (r) , 
with p being the average level density. In fact, it can be shown [51] that the 
expression for the compound reaction cross section 



p _ 1 + k v (2 t + 1 ^ MrM w^ 

*<* ~ k] (2I t + 1)(2/, + 1) ^ + 1} Ea , (Ttf { ) 



can be derived from a Breit-Wigner expression and substituting the resonance 
parameters with average quantities. The Breit-Wigner expression for the cross 
section as the sum of n individual resonances is given by 

p 71-2 i + n r^rf 

^ w = fcf (2/, + !)(]■ + 1) g (2Jfe + 1} (e - Ej+\r^/2y ■ (28) 



The indices % refer to the target, j to the projectile, a to the entrance channel, 
and (3 to the exit channel. The total width of a resonant state k, rj?*, is the sum 
over the partial widths of the individual decay channels rjf* = T£ + rf + . . . 
Using the relation 

/ r^it \ i r r^rf JE1 ^ 2vr rgi* 



^-x^^f(29) 



(£ - £ fc ) 2 + (IT/2) 2 / AE J (E - J4) 2 + (r^/2) 2 A£ r£ 



we can write 



£(2 J + 1)2 W J Z\ J ^ W^U vr) 

J,W \1 J,7T / 

From the last line above it is easily seen that the Hauser-Feshbach cross section 
in Eq. 27 is an averaged Breit-Wigner cross section <t C n = (cbw), when W a/3 = 
1. 

In the general case, the width fluctuation coefficients 
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describe non-statistical correlations between the widths in the channels a and 
/3. They differ from unity close to channel openings and enhance the weaker 
channel [43]. 

Although it might seem tempting to conclude that the cross section of a reac- 
tion proceeding through the compound mechanism should be smooth because 
it is formed from the superposition of amplitudes from a very large number 
of states with random phases, this is a wrong assumption [43]. Ericson [52] 
first showed that the cross sections can continue to show large fluctuations. 
The usual Hauser-Feshbach equations do not account for these fluctuations. 
Therefore, a meaningful comparison to experimental data is only possible af- 
ter averaging the data over a sufficiently wide energy range, comparable to 
the average resonance widths. When using the statistical model to compute 
astrophysical reaction rates (or when deriving rates experimentally directly) 
this is taken care of automatically. However, when using beams with a very 
narrow energy spread it should be noted that the results cannot be directly 
compared to calculations [53]. 

There is a number of available predictions, e.g. [54,55,56,57,58,59,60]. The 
underlying model is the same in all cases [50] but the used nuclear properties 
and prediction of properties are different. Due to the number of reactions to 
be included in reaction networks, global models have to be applied. It should 
be noted that rates derived from the statistical model are not applicable at 
all astrophysical temperatures. The specific temperature limits are given, e.g. 
in [56], and it is advisable to refrain to use those rates below or close to the 
given temperature. 

The most relevant quantities for calculating the statistical model cross sections 
are reaction Q values, optical potentials, level densities, the electromagnetic 
transition strength (mostly El and Ml) for photon transmission coefficients, 
and information on low-lying discrete states. Due to the fact that average 
properties of usually a large number of transitions are used, the results are 
not as sensitive to a variation in the input as the direct reaction cross sections. 
Nevertheless, extrapolations to far from stability can bear considerable errors 
due to effects which are not pronounced close to stability. Modern investiga- 
tions focus on improving the prediction of nuclear properties needed in the 
statistical model. There are basically two types of approaches. First, one tries 
to obtain new data as close as possible to the astrophysically relevant energy 
range for a number of nuclei. From a systematic comparison of predictions 
with various nuclear inputs to the data, theory is tested close to stability. 
Another approach is to implement recent microscopic information to derive 
nuclear properties far off stability and study their impact on the rates. The 
latter approach is the only one possible when dealing with highly unstable 
target nuclei but it is prone to the highest uncertainties. An attempt has to 
be made not just to take nuclear structure predictions at face value but rather 
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Fig. 3. Shown is the nuclear level density ratio of ir = — 1 levels to n = +1 levels 
in 61 Fe as a function of excitation energy. Note that the ratio is considerably lower 
than 1 even several MeV above the neutron separation energy (E scp = 5.58 MeV). 

to learn from comparison of different model and to extract dependences of the 
required nuclear properties. When such basic dependences are known, they 
can be easily implemented, e.g., in a parametrized form, in global reaction 
models. A good example for proceeding this way is the treatment of the nu- 
clear level density. A number of recent microscopic calculations have proved 
that the backshifted Fermi-gas approach [47] perfectly reproduces the shape 
of the level density as a function of energy, both for exotic nuclei and to high 
excitation energies [61,62,63,64,65]. Thus, the determination of the level den- 
sity parameter and the backshift remains. Those can again be extracted in a 
simple manner from globally available calculations of nuclear masses. At low 
excitation energies the Fermi-gas formula has to be replaced by a constant 
temperature formula or a similar approach. At intermediate excitation one 
has to account for the fact that parities are not distributed evenly. Based on 
Monte Carlo shell model calculations a simple model for the partition function 
ratio was suggested recently. We have adopted this model and are currently 
extending it to all nuclei [66,67]. Fig. 3 shows an example for the calculated 
parity ratio in 61 Fe. This will allow not only to establish a better descrip- 
tion of cross sections for nuclei with low separation energies but will also be 
important for checking the applicability range of the statistical model. 

Comparison with recent reaction data has shown that the global statistical 
model approaches work well at stability for reactions involving neutrons and 
protons as projectiles but that there can be problems predicting a capture due 
to the uncertainty in the global optical a potential [47,68,69,70,71,72,73,74,75,76,77,78]. 
The problem lies in the fact that astrophysically relevant ct-energies are close 
to or below the Coulomb barrier and the potentials show a strong energy- 
dependence over this energy range. Thus, the optical potentials are both hard 
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to measure and difficult to calculate. A number of attempts was made to derive 
an improved global description of the potential [53,79,80,81,82,83,84,85,23]. 
Very promising seems the new global potential of [86] which led to good agree- 
ment for a number of reactions [87]. 

A new group of experiments is testing the ground-state photon transmissions 
and thereby the reliability of the used GDR strength [88,89,90,91,92,93,94,95,96,97] 
(see also Mohr et al., this volume). Also there, no discrepancies have been 
found so far. However, concerning the photon strength function it has been 
suggested that there could be some additional strength at low photon energy, 
created by soft-mode vibrations, often called pygmy resonance. A low-lying 
pygmy resonance could increase the neutron capture cross section far off sta- 
bility considerably [17,98]. However, different nuclear structure models still 
disagree about position, strength and existence of this resonance in exotic 
nuclei, e.g. [99,100]. It will only be important if there is additional strength 
below the neutron separation energy, if the statistical model is applicable, and 
if it also appears in nuclei which are sufficiently close to stability so that they 
are not produced in a reaction equilibrium but rather by individual capture 
or photodisintegration reactions. 



4-3 Transition regions with intermediate level densities 

Because of the Coulomb barrier, the Gamow window for charged particle re- 
actions is located at higher energy than for neutron-induced reactions (see 
Sec. 2). Therefore the statistical model of nuclear reactions is applicable for 
the majority of cases with charged projectiles. The only exceptions are proton 
captures close to the proton dripline as encountered in the rp process. The 
situation is different for neutron capture. Although a large number of reac- 
tions at stability proceeds at high level densities, already at stability there 
are some which reach into the regime of intermediate and low level densities. 
As stated above, this happens for nuclei with low level densities at neutron 
shell closures. On the neutron-rich side, unstable nuclei exhibit low neutron 
separation energies soon, leading to formation of the compound nucleus at low 
level density. 

The effective neutron energies are quite low in all astrophysical scenarios, 
E n < 120 keV [68]. Although some scenarios start out from higher plasma 
temperatures leading to higher neutron energies, reaction equilibria are estab- 
lished and individual rates lose importance at such conditions. Since one has 
to cover the energy range of E Q < 120 keV, it becomes obvious that several 
reaction mechanisms might contribute in different parts of that range in nu- 
clei with low level density. Especially the treatment of isolated resonances is 
very difficult because the locations and strengths have to be predicted very 
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accurately in order to apply R-matrix or Breit-Wigner formulas. Furthermore, 
the interplay between competing reaction mechanisms has not been treated 
rigorously so far. It can be shown that there should be no coherent effects 
between direct and compound reactions [II] but there can be interferences 
between single resonances and a direct background [101]. 

The imaginary parts of the optical potentials used might have to be modified 
to account for the flux into other reaction channels. The statistical model uses 
the optical model reaction cross section as given by Eq. (23) to compute the 
formation cross section <7f orm . Thus, it is implicitly assumed that all flux lost 
from the elastic channel is used to form a compound nucleus subsequently 
decaying statistically. In the transition region this is an oversimplification as 
other reaction mechanisms become important, too. Therefore, an increasing 
fraction of the lost flux is accounted for in other reaction channels, not in 
the statistical one, and the imaginary part used for <7f orm would have to be 
decreased accordingly with decreasing level density. A neglection of this effect 
leads to an overestimate of <7f orm and thus of the statistical model cross section. 

Although there have been successful calculations incorporating all reaction 
mechanisms [102,103,104] and explaining differences between resonance and 
activation measurements [105,106] without modification of the potential, this 
conceivable complication should be kept in mind. In consequence, it is very 
important to have measurements directly in the relevant energy range. Extrap- 
olations upwards from lower energies or downwards from higher energies can 
prove treacherous, especially when different reaction mechanisms are involved 
[107]. Even when just measuring elastic data to derive optical potentials, the 
energy range relevant for astrophysics has to be covered. 



5 Conclusion 

The methods to predict nuclear reaction rates are well established. Approaches 
based on nucleon-nucleon interactions have the advantage that they are first- 
principle methods, automatically including all relevant effects such as antisym- 
metrization, correlations, Pauli blocking, etc. On the other hand, depending on 
the model, effective NN interactions might have to be adjusted for each case. 
Common to all microscopic approaches is that the involved numerics limits 
the application to few-nucleon systems. The other class of reaction theories, 
based on optical potentials derived either from microscopic structure consid- 
erations or from phenomenology, average over certain properties and can be 
viewed as approximations to the exact calculation. Their strength lies in the 
wide applicability and the fact that the neglection of certain effects can be 
compensated by slightly adjusted optical potentials (which is automatically 
ascertained when deriving phenomenological potentials from experiment). 
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However, any model can only be as good as the underlying nuclear structure 
model required to provide the necessary input for the reaction calculation. 
While the combination of experimental data and predicted nuclear structure 
leads to good accuracy in the theoretical rates at and close to stability larger 
errors should be expected further off stability. Future radioactive ion beam 
facilities (RIA, GSI upgrade) will provide the possibility for improvement. 
The situation is also eased by the fact that nuclei far off stability are usually 
produced in reaction equilibria (such as the (n,7)-(7,n) equilibrium of the r 
process) where individual reactions are not important. Thus, reactions have 
to be explicitly included only closer to stability (hydrostatic and explosive 
burning, freeze-out from high-temperature scenarios). Due to this fact it is 
actually possible to produce meaningful astrophysical models. We can also 
expect the future facilities to directly study this relevant range of nuclei. 
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